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Abstract — The  ^pography  of  thejwabed  can  strongly  affect  under¬ 
water  sound  pro^gation  in  the  ocean. 'fn-this  rcgartl^aealloor  fea¬ 
tures  fall  into  fheo^overlapping  categories  according  ft  siie:  large 
features  that  block  propagation,  intermediate  features  that  act  pri¬ 
marily  as  sloping  bottoms,  i^d  small-scale  features  that  act  as  scat- 
tertrs.^h^l^is  paper^ statistital  parameters  of  bottom  topography  for 
th^atter  categories ^e  pMsented?  Spatial  wavenumber  spectra 


of  the  seafloor  have  been  estimaied  and  included  in  data 
bases  for  use  in  acoustic  modeling.  Examples  of  this  are  the 
“geoacoustic  models"  of  Hamilton  (8)  based  on  I)  w  situ  and 
.  laboratory  measurements  on  sediments.  2)  seismic  experiments, 
and  3)  acoustic  experiments. 

The  effect  of  seafloor  topography  on  underwater  sound 


of  ocean  bottom  and  subbottom  roughness  are  determined  from  oar-  f  propagation  is  a  function  of  experimental  geometry  and 
row-beamwidth  echosounding  and  seismic  reflectiM  er^g.  The./  frequency.  Topographic  features  fall  into  three  overlapping 
spectra  are  compared  to  the  expression  /’(Jf  I  =  CK  "*1,  where  PiKt  is  -  .  , ,  i  r  , 

.r  .  1 -i  1,  i-.  .  ,  .4.  categories  1) large  features  that  block  propagation, 

the  power  spectral  density,  C  is  a  proportionality  constant,  A,  is  the  /  e  '■  r 

wavenumber,  and  his  a  constant  that  characteriies  the  clau  of  rough-  2)  intermediate  sized  features  that  primarily  act  as  sloping 
ness.  The  parameter  b  is  often  assumed  to  be  3;  however, <be  preoewt  ■  bottoms,  and  3)  small-scale  features  that  act  as  scatterers. 


if  study  shows  that  ft  can  range  from  about  I  to  S.  Topographic  samples 
were  found  to  have  probability  density  functions  which  were  both 
non-(iaussian  and  Gaussian.  It  is  suggested  that  a  first^rder  rough¬ 
ness  data  base  include  band-limited  mean  square  tVMSirougb- 
ncss;  A' j  and  Jf  Vfthe  wavenumbers  of  the  estimatel;  ft;  sediment  type; 
physiographic  province,  water  depth,  and  location.  ;  ‘c.  "  f 

-  INTRODUCTION 

SEAFLOOR  ROUGHNESS  is  an  important  factor  in  acous¬ 
tic  propagation.  Properties  of  roughness  are  not  only  a 
means  for  studying  seafloor  geology  (14j,  (161,  but  also 
provide  a  method  for  seafloor  classification  (S|.  This  paper 
deals  with  a  quantitative  description  of  seafloor  topography 
for  use  in  acoustic  problems.  First,  various  types  of  roughness 
parameters  that  have  been  used  as  input  to  acoustical  models 
are  reviewed.  Then  estimates  of  seafloor  and  subbottom 
roughness  obtained  from  stabilized  narrow-beam  echosound- 
ings  arc  presented.  These  data,  along  with  data  presented 
from  the  literature,  can  prosdde  interim  estimates  of  roughness 
parameters  until  an  extensive  roughness  data  base  is  estab¬ 
lished.  Finally,  the  form  of  a  first-order  seafloor  rougliness 
data  base  is  suggested. 

SEAFLOOR  TOPOGRAPHY  AND  UNDERWATER 
ACOUSTICS 

The  interaction  of  sound  with  the  seafloor  depends  upon 
botioni  dciisii) .  sound  attenuation,  sound  velocity,  and  inter¬ 
face  rougliness.  The  densii)  .  sound  velocity,  and  attenuation 
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Only  topography  of  the  first,  and  to  some  extent  the  second 
categories,  is  readily  available  for  use  in  acoustic  modeling. 

Deterministic  Data  Bases 

Topographic  features  of  the  seafloor  of  the  first  and  second 
categories  given  above  may  be  described  deterministically 
and  input  iiiio  range-dependent  acoustic  propagation  models 
such  as  Parabolic  Equation  |22]  and  CRASS  |4]  .Topographic 
data  are  usually  obtained  from  bathymetric  charts  or  data 
bases  such  as  SYNBABS.  a  computerized  bathymetric  data 
base  and  software  system  that  synthesizes  great-circle  bathy¬ 
metric  profiles  from  average  depth  in  1/12  degree  cells  |23l . 

Recently  developed  single-interaction  scattering  models, 
sucli  as  Facet  Ensemble  (1131 .  (18])  require  input  of  a  high- 
rcsoluiion  topographic  profile.  Such  profiles  are  difficult  to 
obtain  on  a  global  scale,  but  a  data  base  to  support  such 
modeling  might  consist  of  a  scries  of  profiles  from  areas  in 
which  there  is  unifonn  small-scale  roughness  (roughness 
provinces)  The  data  base  could  be  either  a  part  of,  or  separate 
from,  the  statistical  data  bases  described  below. 

Statistical  Data  Bases 

Intermedialc-scale  and  small-scale  features  cause  scattering 
of  sound  and  errors  in  range  and  bearing  estimates  (11]. 
Different  statistical  parameters  of  roughness  are  required 
for  different  scattering  theories.  Eckart  [6]  has  shown  the 
spatial  wavenumber  spectrum  to  be  an  important  factor  in 
the  scattering  of  sound  from  a  randomly  rough  surface. 
Clay  ct  al  |3|  showed  that  the  cuherent  component  of  the 
specularly  scattered  sound  is  sensitive  to  the  probability  density 
fiinctioii  (I'lJl  )  of  the  displacements  of  the  rough  surface. 
For  the  case  ol  a  (uussian  PDF,  F.ckart  |6|  showed  that  the  co¬ 
herent  coiiiponciit  of  sound  reduces  to  a  simple  expression  in - 
viilvmg  the  Root  Mean  Square  (R.MSi  loiighness  of  the  surface. 
I  here  have  been  suggestions  that  the  seafloor  roughness 
PDJ 's  tend  to  he  approximately  Gaussian  1 15] . 

In  order  to  represent  sfiatial  wavenumber  spectra  for  areas 
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of  the  seafloor,  a  simplified  model  is  convenient.  Spectra  may 
be  approximated  by  the  expression =  CA^”*,  where /XA') 
is  the  spatial  wavenumber  power  spectral  density,  C  is  a 
porportionality  constant,  K  is  the  spatial  wavenumber,  and  b  is 
a  constant  that  is  characteristic  of  the  class  of  roughness 
(analogous  to  noise  class,  e.g.,  white  noise.  Brownian  noise). 
Nye  {19)  has  used  a  dimensional  analysis  to  demonstrate 
that  for  the  case  of  b  =  3,  the  units  of  spatial  wavenumber 
power  spectral  density  (meters  cubed)  cancel  and  the  topog¬ 
raphy  appears  to  have  the  same  roughness  for  all  scales.  An 
example  of  the  acoustic  significance  of  b  is  shown  by  Marsh’s 
[17]  theory  of  scattering  from  a  totally  reflecting  randomly 
rough  surface.  For  Fixed  grazing  angle,  the  backscattering 
coefficient  varies  as  where  k  is  the  acoustic  wavenum¬ 

ber.  Note  that  for  b  =  3,  the  backscattering  would  be  inde¬ 
pendent  of  acoustic  frequency. 

Most  studies  of  seafloor  topography  have  been  qualitative. 
A  few  quantitative  studies  have  dealt  witli  very  large-scale 
topography.  For  acoustic  analysis,  statistical  parameters 
relating  to  the  roughness  of  the  acoustic  interaction  zone 
are  required.  In  the  next  sections,  the  roughness  important 
to  acoustic  interaction  is  discussed. 

SEAFLOOR  ROUGHNESS  APPLICABLE  TO  LOW-  AND 
MEDIUM-  FREQUENCY  SOUND 

The  area  of  interaction  for  sound  reflecting  from  the 
seafloor  may  be  estimated  by  the  size  of  a  Fresnel  Zone.  For 
lOQ-Hz  sound  (acoustic  wavelength  =IS  m)  and  a  20“  grazing 
angle  and  with  surface  source  and  surface  receiver  in  a  4000-m 
ocean,  the  dimensions  of  the  First  Fresnel  zone  calculated  by 
Kerr’s  {I2j  method,  are  1900  m  by  600  m.  By  the  Rayleigh 
criterion  {21),  the  heights  of  roughness  within  the  First 
Fresnel  Zone  must  be  greater  than  about  A/(8  sin  0),  where  X 
B  the  acoustic  wavelength  and  0  is  the  grazing  angle,  or  S.5  m 
to  appear  as  a  “rough"  surface  to  incident  sound.  To  delineate 
seafloor  features  of  this  scale  requires  better  resolution  than 
conventional  wide-beam  echosounders  can  offer;  they  com¬ 
monly  have  a  60*  beamwkith,  which  would  imply  a  4600-m 
diameter  ensoniFied  area  for  an  ocean  depth  of  4000  m. 

One  method  of  achieving  the  required  solution  is  to  use 
a  stabilized,  very  narrow-beam  echosounder  [7] .  Data  from 
this  type  of  echosounder  were  obtained  by  using  the  beam 
of  highest  resolution  of  the  stabilized  12-kHz  multibeam 
array  sonar.  Depths  obtained  from  the  center  beam  (normal 
incideiKe)  were  determined  to  1  m  by  precise  measurement 
of  the  sound  traivl  time.  The  ensonifled  area  (to  the  -3  dB 
point)  was  less  than  90  m  in  diameter  and  adjacent  samples 
did  not  overlap  because  the  sampling  interval  was  about  100 
m.  Sample  series  of  topographic  data  were  adjusted  to  »  zero 
mean,  and  passed  through  a  high  pass  spatial  Filter  (low  cut 
wavenumber  0.003  m~').  Probability  density  functions  and 
power  density  spectra  were  then  computed  from  the  Filtered 
data.  A  2048-point  diKrete  Fourier  transform  with  a  Hann 
Window  was  used  to  obtain  raw  spectral  estimates.  Averages 
of  10  adjacent  estimates  were  used  to  produce  a  smoothed 
spectral  estimate  having  a  resolution  of  0.0003  m"'  for  the 
band  up  to  0.03  m~'. 

To  reduce  the  effects  of  system  noise,  navigational  un¬ 


certainties,  and  heave,  only  data  obtained  under  optimum 
conditions  are  used  for  this  study.  Aliasing  may  affect  a  spec¬ 
trum  if  substantial  energy  occurs  at  frequencies  higlier  than 
tlie  spatial  sampling  frequency  (I  sample  per  100  m).  How¬ 
ever.  as  the  beamwidth  of  the  echosounder  is  not  infinites¬ 
imally  small,  the  measurement  system  may  act  as  an  antialias¬ 
ing  Filter.  Other  processing  effects  include  bias  due  to  leakage 
from  one  band  to  another.  Leakage  effects  have  been  mini¬ 
mized  by  using  appropriate  windows. 

By  using  the  same  measurement  system,  processing,  and 
estimation  techniques  on  a  wide  range  of  seafloor  types  (Table 
I).  First-order  estimates  of  the  probability  density  functions 
(Fig.  1)  and  power  density  spectra  (Fig.  2)  can  be  obtained. 
These  functions  and  the  RMS  roughness  arc  band-limited 
parameters,  since  they  pertain  to  the  band  of  topographic 
wavenumbers  sampled  by  the  measuring  system  and  the  high- 
pass  processing  Filter.  The  PDF’s  appear  to  have  both  Gaussian 
and  non-Gaussian  distribution.  The  values  of  h  (Table  I) 
were  obtained  by  a  logarithmic  least-square  Fit  of  each  power 
spectrum  for  those  values  that  were  above  measurement  sys¬ 
tem  noise  (Fig.  2).  These  b  values,  which  vary  from  about  1 
to  S,  have  a  greater  variation  than  that  reported  by  other 
investigators.  Nye  (19)  concluded  that  spatial  wavenumber 
spectra  of  widely  different  types  of  land  topography  have 
the  approximate  form  corresponding  to  6  =  3,  even  though 
the  values  of  C  are  greatly  different.  Marsh  |17j  compiled 
power  spectra  of  nine  topographic  surfaces,  including  four 
sea  bottoms  which  followed  the  form  corresponding  to  6  =  3 . 
The  lake-bottom  spectrum  reported  by  Horton  ei  al.  (10) 
has  the  form  6  =  0.  Bell  (1)  calculated  spectra  for  North 
Pacific  abyssal  hills  from  numerous  sources  including  deep- 
tow  echosounding  data  and  found  that  b  varies  from  2.0 
to  2.S  for  wavelengths  less  than  10  km  and  that  b  was  about 
1.0  for  longer  wavelengths.  The  wide  range  in  b  for  the  sea¬ 
floor  is  not  unexpected,  since  there  are  many  unrelated 
processes  that  act  to  form  the  relief.  This  is  in  contrast  to  the 
constant  value  of  6  =  3  for  the  equilibrium  range  of  wave- 
number  spectra  of  fully  developed  wind-blown  sea  surfaces 
(20),  where  roughness  results  from  a  single  mechansim.  One 
characteristic  of  all  seafloor  power  spectra  for  vitually  all  scales 
of  topography  is  that  b  is  rarely  less  than  I ,  indicating  that 
the  power  is  concentrated  in  the  longer  wavelengths.  This 
suggests  that  features  that  are  tall  relative  to  their  horizontal 
dimensions  are  rare.  Such  features  would  tend  to  be  unstable 
and  short-lived  in  the  ocean  environment. 

RMS  roughness  estimates  determined  in  this  study  (Table 
i)  are  consistent  with  those  found  by  Clay  and  Leong  (2). 
Further,  physiographic  provinces  appear  to  be  characterized 
by  certain  ranges  of  RMS  roughness.  However,  there  is  no 
apparent  relationship  in  these  data  between  b  and  RMS 
roughness  or  b  and  physiographic  province.  It  appears  that 
additional  studies  with  much  larger,  higher  resolution  data  sets 
are  required  to  determine  if  there  are  relationships  between 
b  and  seafloor  type. 

If  roughness  spectra  can  be  approximated  by  the  exponen¬ 
tial  expression,  a  statistical  data  base  might  include  parameters 
such  as  band-limited  RMS  roughness;  K,  and  Kj  wavenumber 
bounds  of  the  estimate;  6;  sediment  type;  water  depth:  physio- 


> 


50 


IEEE  JOURNAL  OF  OCEANIC  ENGINEERING,  VOL.  OE-9.  NO.  I.  JANUARY  1984 


TABLE  I 

RMS  ROUGHNESS  AND  SPECTRAL  SLOPE  PARAMETER  6  OF 
BANO-LlMITEO  TOPOGRAPHY 
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Spatial  wavenumbeis  arc  0.003  to  0.03  m~l  except  those  denoted 
by  *,  where  spatial  wavenumbers  are  0.00006  to  0.003  m-i.  In  cases 
where  the  roughness  is  less  than  the  resolution  of  the  measuring  system, 
the  upper  limit  of  RMS  roughness  is  given  and  b  is  not  estimated. 


Fig.  I.  Probability  density  function  of  ocean  bottom  topography.  The 
solid  lines  represent  the  Gaussian  PUi  ol  topography  filtered  by  a 
high-pass  spatial  filter,  wavcnumbei  0.003  m  l.  (a)  Norwegian 
Sea-  Basin,  (b)  Norwegian  Sea-Marginal  Plateau. 
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I'ig  3.  Spatial  wavenumhei  power  cpeeiral  density  ol  sejlloor  topog¬ 
raphy.  The  solid  lines  represent  the  logarithiiiK'  regresMun  for  spec¬ 
tral  values  that  arc  above  the  measureiiient  noise  .(a)  Norwegian 
Sea  Basin,  lb)  Norwegian  Sea  Marginal  Plateau 

yraph  ptrivincc.  and  gcrigrapliic  location  Tlic  KMS  roughness 
lor  other  Ireqiiciicy/waveiiiiniber  bands  could  be  estimated 
from  the  exponential  expression,  l  uilhct  work  is  required 
to  determine  if  such  a  relatively  simple  data  base  can  ade¬ 
quately  represent  boiUHit  roughness. 
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At  vety  low  fiequencies  (less  than  20  11/ 1.  acoustic  wave¬ 
lengths  arc  long  (greater  than  75  m>.  ihe  elicctivc  attenuation 
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low,  and  the  Fresnel  zone  size  large  (dimensions  proportional 
to/^'^*).  A  significant  amount  of  very  low-frequency  energy 
can  pass  through  the  water-sediment  interface  and  interact 
with  the  subbottom.  Whether  the  water-sediment  interface 
or  a  subbottom  interface  is  the  principal  scattering  surface 
will  depend  on  experimental  geometry,  roughness  of  the 
interfaces,  acoustic  wavelength,  sediment  thickness,  sediment 
density,  and  the  sound  attenuation  and  sound  velocity  in  the 
sediment.  For  large  areas  of  the  world's  oceans,  the  principal 
subbottom  interface  for  VLF  sound  is  the  sediment-basalt 
interface. 

To  obtain  statistical  properties  of  both  the  water-sediment 
and  the  sediment-basalt  interfaces,  large-scale  roughness 
data  were  obtained  from  seismic  reflection  records  and  wide- 
beam  echosounding  data  from  the  North  Atlantic  and  North 
Pacific  Oceans.  The  seismic  reflection  records  provide  data 
for  both  interfaces.  The  echosounding  records  were  made  in 
areas  of  little  or  no  sediment  cover.  Depths  were  determined 
at  intervals  of  I  km  along  profiles.  The  accuracy  of  these 
data  varied  between  S  and  30  m,  depending  on  the  seismic 
recording  configuration.  The  seismic  interface  depths  were 
then  adjusted  for  a  constant  sediment  velocity  layer  (1.6 
km/s).  While  these  data  do  not  have  the  resolution  of  the 
narri'w-bandwidth  data  mentioned  in  the  previous  section, 
they  provide  an  estimate  of  rougliness  in  the  0.00006  to  0.(X)3 
m~ '  wavenumber  band.  The  roughness  of  the  smaller  wave- 
numbers  of  this  band  is  applicable  to  scattering  at  the  lower 
frequencies  of  VLF  sound. 

Power  density  spectra  and  probability  density  functions 
were  obtained  from  the  digitized  data  that  have  been  de¬ 
trended  and  have  had  the  mean  removed  Table  1  shows 
estimated  large-scale  roughness  (0.(X)006  to  0.003  m~' 
wavenumbers)  statistics  for  the  sediment-basalt  interface  of 
the  North  Atlantic  and  North  Pacific  Oceans.  These  estimates, 
based  upon  50  sample  profiles  from  each  ocean,  yielded  mean 
values  of  1.8  and  1.6,  respectively,  for  parameter  b.  This  is 
consistent  with  results  reported  by  Bell  |l],  who  found  that 
b  for  ihe  slope  of  the  North  Pacific  abyssal  hills  was  2.0  to 
2.5  for  wavelengths  less  than  about  40  km  with  a  lower  value 
for  longer  wavelengths.  Bell’s  results  showed  a  large  apparent 
scatter  with  only  a  few  data  points  in  the  long  wavelength 
range.  The  standard  deviation  resulting  from  fitting  the  ex¬ 
ponential  approximation  to  each  power  spectrum  was  0.4, 
which  iscompaiablc  to  that  found  by  Beil  |l| . 

The  PDF  of  all  samples  free  of  seamounts  and  fracture 
zones  were  found  to  have  a  distinct,  generally  symmetric, 
central  tendency.  This  is  in  agreement  with  the  findings 
of  Krause  ct  al  |15)  for  the  North  Pacific  and  Holcombe 
19|  for  the  North  Atlantic. 

While  the  values  of  b  for  basalt  are  similar  for  both  oceans, 
the  average  RMS  roughness  is  significantly  different  (Table  1). 
RMS  roughness  for  the  basaltic  basement  of  the  two  oceans 
in  the  spatial  wavelength  range  from  5  to  100  km  is  estimated 
to  be  259  ♦  74  m  for  the  North  Atlantic  and  99  ±  36  ni  for 
the  North  Pacific  Oceans.  This  analysis  excluded  the  large 
fraciuic  /ones.  The  uncertainty  indicated  is  one  standard 
deviation  of  the  individual  estimates.  Ihe  means  are  distinctly 
diflcient,  with  the  North  Atlantic  having  the  higher  value. 
Holcombe  |9|  has  estimated  mean  rebel  m  the  North  Atlantic 
by  hand  tabulation  of  peak-to-vallej  heights  and  by  averaging. 
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Fi*.  3.  Total  sediment  thickness  veisus  ratio  of  chanite  in  RMS  rouith- 
ness  to  basenieni  RMS  loufihness. 


An  approximation  comparison  of  these  two  results  can  be 
made  by  assuming  the  topography  to  be  sinusoidal.  Hits  as¬ 
sumption  yields  a  conversion  of  the  Holcombe  peak  to  peak 
amplitude  to  RMS  roughness  of  244  i  60  m.  Tlie  value  is 
in  agreement  with  our  data.  Increasing  sediment  cover  de¬ 
creases  relief  at  the  water-sediment  interface.  As  shown  in 
Fig.  3.  the  reduction  in  RMS  roughness  is  approximately 
proportional  to  sediment  thickness;  however,  there  is  con¬ 
siderable  scatter  in  the  data. 

CONCLUSION 

Statistical  properties  of  seafloor  and  subboitom  interfaces 
have  been  calculated  for  a  variety  of  seafloor  types  and  loca¬ 
tions.  Conclusions  of  this  study  are.  l)the  roughness  slope 
parameter  b  varies  from  about  I  to  S:  2)  while  many  proba¬ 
bility  density  functions  approximate  a  Gaussian  distribution, 
there  are  exceptiofu;  3)  the  difference  between  the  RMS 
roughness  of  the  basaltic  basement  and  the  overlying  sedi¬ 
ments  is  approximately  proportional  to  sediment  thickness; 
and  4)  the  RMS  roughness  of  the  basaltic  basement  is  much 
larger  in  the  North  Atlantic  than  the  North  PaciHc. 

With  the  increased  availability  of  stabilized  multibeam 
echosounders,  there  is  a  potential  for  developing  large  data 
bases  of  high-resolution  seafloor  topography  sutistics.  We 
have  suggested  that  a  (Ust-order  roughness  dau  base  include 
the  following:  band-limited  RMS.  1C|,  K^.b,  sediment  type, 
physographic  province,  water  depth,  and  location.  An  addi¬ 
tional  data  base  might  also  include  actual  topographic  samples 
representative  of  the  various  seafloor  types.' 
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